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ABSTRACT 

This paper examines the allowed amount of IG (intergalactic) dust, which 
is constrained by extinction and reddening of distant SNe la and thermal 
history of IGM (intergalactic medium) affected by dust photoelectric heating. 
Based on the observational cosmic star formation history, we find an upper 
bound of x, the mass ratio of the IG dust to the total metal in the Universe, as 
X < 0.1 for 10A < a < O.lum and X < 0.1(a/0.1 \aa) for O.l^m < a < lum, 
where a is a characteristic grain size of the IG dust. This upper bound of 
X ~ 0.1 suggests that the dust-to- metal ratio in the IGM is smaller than 
the current Galactic value. The corresponding allowed density of the IG dust 
increases from ~ 10~ 34 g cm" 3 at z — to ~ 10~ 33 g cm" 3 at z ~ 1, and 
keeps almost the value toward higher redshift. This causes IG extinction of 
< 0.2 mag at the observer's B-band for z ~ 1 sources and that of < 1 mag for 
higher redshift sources. Furthermore, if E{B — V) ^0.1 mag at the observer's 
frame against z>l sources is detected, we can conclude that a typical size of 
the IG dust is < 100 A. The signature of the 2175 A feature of small graphite 
may be found as a local minimum at z ~ 2.5 in a plot of the observed E(B~V) 
as a function of the source redshift. Finally, the IGM mean temperature at 
z < 1 can be still higher than 10 4 K, provided the size of the IG dust is < 100 

A. 
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1 INTRODUCTION 

As long as there is dust between radiation sources and observers, the dust extinction and 
reddening 1 must be corrected if we want to realize the nature of the sources. We should 
examine how much extinction a nd reddening there ar e. The extinction property in the Galaxy 



is a well studied example (e.g. 



Draine fc Lee 



1984). Using HI and far-infrared emission as 



tracers for the dust column density, we can obtain the extinction amount by the Galactic 



dust with reasonable accuracy (jBurstein fc HeileiJl982l : ISchlegel. Finkbeiner. fc Davislll998h . 



Although the dust distribution and properties in the external galaxies are not well known 
yet, we can correct the dust extinctio n in the galaxies by using some indicators of the 
extinction with some assumptions (e.g.. 
intergalactic (IG) dust? 



Buat et al 



2002) . How about the extinction by the 



We have already known the fact 
even at redshift larger than 3 (e.g., 



;hat some metal e l ements exist in th e Lyman a clouds 



Cowie et al 



1995 



Telfer et al 



200J). It suggests that 



the dust grains also exist in the low-density intergalactic medium (IGM). Such diffuse 
IG dust causes the IG extinction and reddening, which may affect on our understanding 



of the Universe significantly. One might think that the IG dust amount is negligib 



cause such a significant IG reddening is not reporte c 



Gheng. Gaskell. fc Koratkar 



1991 



Riess et al 



1998; 



in the previous studie s rtTakase 



Perlmutter et al 



e be- 



1999). However, the 



wavelen g th de pendence of the IG extinction is quite uncertain. If it is gray as suggested by 
Aguirrd (J1999J), a large extinction is possible with no reddening. Nobody can conclude that 
the IG dust is negligible because of no observable reddening. 

Theoretically, it is predicted that metals sy nthesized in supernova (SN) explosions form 



2001 



Nozawa et al. 



2003 



Schneider. Ferrara. fc Salvaterr. 



into the dust grains in the cooling ejecta of SNe dKozasa fc Hasegawall987 



Todini fc Ferrara 



20031 ). Recently, t hermal emis- 



sions of such dust from two supernova remnants, Cas A and Kepler, are detected ((Dunne et al 



2003; 



Morgan et al 



2003). In a very high-z universe, SNe of massive Population III stars 



formed in low mass halos, which are likely to b e the main site of the first star form ation, 
can disperse the produced metals into the IGM (|Bromm. Yoshida. fc Hernquistll2003h . The 
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1 In this paper, we call the total absolute amount of the absorption and scattering at a wavelength just extinction, and the 
differential extinction between two wavelengths reddening. 



dust grains may be also d ispersed into the 



even in a z > 10 universe (jElfgren fc Desert 
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GM. Therefore, the IG dust grains may exist 



2003). 



Extinction by the IG dust may affect on the determination of the cosmological parame- 
ters from observations of SNe. The observed dimming beyond the geometrical dimming in 



the empty 



space of distant 



constant ( Riess et al 



1998; 



0.5) Type la SNe , which is attributed to the cosmological 



Perlmutter et al. 



tinction without the cosmological constant (j Aguirre 



1999), can be reproduced by the gray IG ex- 



1999). iGoobar. Bergstom. & Mortse 



Riess et al 



2Q0J 



( 20021 ) show that the apparent brightening of the farthest SNe la (z = 1.7 
can be also explained by the gray IG extinction with zero cosmological constant if the dust- 
to-gas ratio in the IGM decreases properly with increasing redshift. Therefore, to know the 
amount of the IG dust is also important in the cosmological context. 

The evidence of the IG dust should be imprinted in the cosmic microwave background 
(CMB) and infrared background because the dust emits thermal radiation in the wave 



band from the far-infrared to submillimetre 



Wright 



1981 



'submm) (jRowan-Robinson. Negroponte. fc Silk 



Elfgren fc Desertl 120031 ) . Although the COBE data provides us with 



only a rough upper bo und on the IG dust 



Aguirre fc Haimanl 120001) . the data of WMAP (jSnergel et al 



Loeb fc Haima: 



1997 



Ferrara et al. 



1999: 



20031 ) may be promising. The 



subm m background radiation will give a more strict constraint on the IG dust (j Aguirre fc Haiman 



2000). 



Recently, we hav e proposed a new constr aint on the IG dust amount by using thermal 



history of the IGM 



Inoue fc Kamava 2003). Si nce t he dust photoelectric heating is very 



efficient in the IGM (N ath. Sethi, fc Shchekinov 



19991 ) . the theoretical thermal evolution of 



the IGM taking into account of the heating by dust violate s the observationa l temp erature 



evolution if too much IG dust is input in the model. Hence, llnoue fc Kamaval (j2003h obtain 
an upper bound of the IG dust amount in order that the theoretical IGM temperature should 
be consistent with the observed one. The obtained upper bound of the dust-to-gas ratio in 
the IGM is 1% and 0.1% of the Galactic one depending on the IG grain size of ~ 100 A- 
0.1 /im and ~ 10 A, respectively, at redshift of ~ 3. 

In this paper, with help of distant SNe la observation, we extend our previous approach 
in order to discuss the upper bounds of the IG dust extinction and reddening. In the next 
section, we start from the cosmic star formation history (SFH) to specify the IG dust amount 
at each redshift. According to the assumed SFH, we can estimate IG dust extinction and red- 
dening at each redshift theoretically. In section 3, we comment on observational constraints 
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from the extinction and reddening of distant SNe la. In section 4, further constraints are 
presented by comparing theoretical and observational thermal histories of the IGM. Based 
on the allowed amount of the IG dust, we also discuss some implications from our results in 
section 5. The achieved conclusions are summarised in the final section. 

Throughout the paper, we stand on a A-cosmology. That is, we constrain the amount of 
the IG dust in order that the IG dust should not affect on the determination of the cosmo- 
logical parameters fro m distant SNe la. This is because the flat universe is fav ored by results 



of CMB observations ( Jaffe et al 



2001; 



Prvke et al 



the matter cannot make the flat universe (Percival et al 



20021 



Spereel et al.ll2003T) . whereas only 



20011 1 . Furthermore, the recent 



observations of the X-ray scattering halo around high-z QSOs suggest too small amount 



of the IG dust to expl ain all amou nt of the dimming of the distant SNe la (jPaerels et al 



2002; 



Telis et al 



2002, but see also 



Windtl 



20021 ). iMortsell Goobarl ifenoah also reach the 



same conclusion by analyzing the observed colours of the SDSS (Sloan Digital Sky Sur- 
vey) quasars. The following cosmological parameters are adopted: H = 70 km s _1 Mpc -1 , 
Q M = 0.3, Q A = 0.7, and Q h = 0.04. 



2 STAR FORMATION HISTORY AND INTERGALACTIC DUST 

To estimate the IG extinction and reddening theoretically, we must investigate production 
of dust at each redshift. Since dust is made of metals, a cosmological evolution of the metal 
amount should be specified. As metals are products of stellar evolution, therefore, we shall 
specify t he cosmic SFH as a first step. 



Since 



Madau et al. 



(1996|), researches of the cosmic SFH are extensively performed. In 
figure 1, we show observational star formation rates in a unit comoving density as a function 
of redshift. The cross and open symbols are estimated from the Ha line and the rest-frame 
ultra-violet (UV) luminosities not corrected by the interstellar dust extinction. Hence these 
are lower limits. The filled-circles are estimated from the submm data. Due to the small 
statistics, the uncertainty of the submm data is rather large. The real SFH is still uncertain 
because we do not know the suitable correction factor against the internal dust extinction. 
In this paper, therefore, we adopt two example models: high and low SFHs, which are shown 
in figure 1 as solid and dashed curves, respectively. For the low SFH case, we have employed 
a conservative correction factor for the internal dust extinction. The high SFH case is set to 
be compatible with the submm data. Quantitatively, these SFHs are formulated as 
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Figure 1. Cosmic star f ormation rate densit y as a function of r edshift. The plotted po ints with error-bars are estimated 
from the Ha data: cross iGallego et alJll995l) and inc lined cross if Tresse fc Madded ll99Sl) : fr om the rest-frame ultra-violet 
data: open-triangles iLillv et alJl996ft open-diamondsJConnollve^al]^997l) , open-pentagons iMad au, Pozzetti, & Dickinson 

circles iSteidel et aJjl99Sft . open-squares llffiavaliscoet^Ll20o3) . open-in verse-pentagon jfeat^^™ . _ 
inverse-triangle iBouwens et all2003l) : and from the submillimetre data: filled-circles fearger. Cowie. fc Richards 2000). These 
points are adjusted to the cosmological parameters assumed in this paper. The solid and dashed curves are the models of high 
and low star formation histories, respectively. 



1998), open 



2003), open- 



Psfr(^) 



3.3 



i±2 ' (for z<:i) 



O.1M yr" 1 Mpc 



-3 



2 



-1.5 



(for z > 1 low SFH) 
(for z > 1 high SFH) 



(1) 



where Psfr i s the star formation rate per unit comoving den sity. Recently, some authors 



suggest that a SFH like the high case is more likely (e.g., ISpringel &: Hernquist 



Giavalisco et al. 



2003; 



20041 ). Hence, we discuss the high SFH case mainly in the following. 



Once a SFH is specified and if the instantaneous recycling approximation is adopted, the 
cosmic mean metallicity 2 evolution is determined by: 



Z(z) 



Vz 



2S 



dz' 



n / Psfr( z ) tt( f \n , a ' (^) 

il h p cfi Jz H{z'){l + z') 

where yz is the produced metal mass fraction for a unit star forming mass, p Cj o is the 
current critical density, zo, is the starting redshift of the cosmic star formation, and H{z) is 
the Hubble constant at the redshift z. We have assumed yz to be a constant for simplicity. 



If the Salpeter initial mass function (0.1-125 M ) is assumed, yz = 0.024 ()Madau et al 



19961 ) . We assume z$ = 10 in this paper. This does not affect on the results obtained in 
the following sections because the measure of time along the redshift is small in the high-;? 



2 In this paper, wc define metallicity as the mass fraction of elements heavier than Li as the usual way. 
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Figure 2. Mean cosmic metallicity evolution. The vertical axis is normalized by the Solar metallicity Zq = 0.02. The solid 
and dashed curves are the cases of the high and low star formation histories, respe ctively. The s tartin g redshift of the cosmic 
star formation is se t to be zg = 10. The observed range of oxyg en abun d ance by ITelfer et alJ l|2p023, carbon abundance by 
ISchave et alj l2003t) . and lower limit of metallicity obtained by lSongailal 1200 ll) are overlaid as the shaded squares and the 
upward arrow, respectively. We note that lAguirre et al 1 J2004) find [Si/C]~0.8. 

universe. Indeed, the cosmic metallicities in z < 3 for various z$ become nearly equal to 
each other if z s > 5. 

In figure 2, we show the mean cosmic metallicity as a function of redshift. In t he figure. 



the lo wer limit of metallicity in the Lya clouds measured 
()200lt ) and the ranges of carbon ab u ndanc e obtained by 



abundance obtained by 



Telfer et al 



rom C IV an d Si IV by 



Schave et al. 



Songaila 



(2003J) and oxygen 



( 20021 ) are overlaid. Our theoretical estimate of total 
metallicity in the universe is compatible with the observed oxygen abundance in the IGM. 
This may mean that a large part (~ 50%) of metal produced in galaxies exist out of galaxies. 
On the other hand, our estimate is much larger than the observed carbon abundance. This 
may indicate that carbon is not suitable to trace the cosmic mean m etallicity in the early 
phase since the first metal pollution is made by the Type II SNe (e.g.. 



Page! 



1997). Indeed, 



Aguirre et all j20pj ) find [Si/C]~ 0.8. 

Let us introduce one parameter to describe the amount of the IG dust; the ratio of the 
IG dust mass to the total metal mass in the universe defined as 
IG dust mass V lGM 



X 



(3) 



total metal mass Z 
where £) IGM is the dust-to-gas mass ratio in the IGM. In principle, this parameter is de- 
termined by the transfer mechanism of dust grains from galaxies into the IGM. Although 



some authors have tried to solve this problem approximately (e.g.. 



Ferrara et al. 



1990 



1991 
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Aguirre et al.ll200lh . the results are not still conclusive. This is because we must solve prob- 



lems of the magneto-radiation-hydrodynamics of dusty plasma finally. Here we approach the 
parameter \ by another way; we constrain the parameter by using the observational data of 
distant SNe la (§3) and the thermal history of the IGM (§4). While the parameter x may 
evolve along redshift, we treat it as a constant for simplicity. Hence, the obtained values of 
X in the following sections are regarded as those averaged over redshift. 

A dust model must also be specified. In this paper, we adopt the "graphite" a nd th e 



"smoothed astronomica 



silica te" models by 



Draine fc Leel (|l984l) : lLaor fc Drainel (|l993l) : 



Weingartner fc Draind ([2001al ). which can explain the interstellar dust properties in the 



Galaxy and the Magellanic clouds very well. Although there is no evidence that the IG dust 
is the same as the Galactic dust, we assume them as a working hypothesis. 

The grain size distribution in the IGM is also quite uncertain. As a first step, we assume 
all grains to have the same size. This means that the grain size in the current paper indicates 
a characteristic size of the IG dust (i.e. an averaged size by a certain way). lAguirrel ((1999) 
suggests a selection rule in the transfer of dust grains from the host galaxies into the IGM; the 
small grains are destroyed by the thermal sputtering process when the grains are transfered 
through the hot gas halo of the host galaxies. Indeed, this theoretical suggestion is very 
interesting to realize the gray dust model. However, we also examine the possibility of the 
small IG grain because the suggestion of the selection rule is not confirmed observationally 
at the moment. 



3 CONSTRAINT FROM OBSERVATIONAL DIMMING OF 
SUPERNOVAE 

In this section, we constrain the amo unt of the IG dust by means of observational dimming 



of distant Type la SNe. According to 



Riess et al 



(200l|), the dimming of SNe la at z ~ 0.5 



is ~ 0.2 ± 0.1 mag against the empty universe. Our policy is that the IG dust extinction 
does not affect the interpretation of the cosmological constant. Thus, we attribute the 0.2 
mag dimming of SNe la at z ~ 0.5 to the cosmological constant and regard the quoted error 
(0.1 mag) as an upper limit of the contribution by the IG dust to the dimming. Although 
the cosmological dimming does not depend on the observed wavelength, the dimming by 
the IG dust may depend on the wavelength. As the distant SNe la are observed in B and 
y-bands and 5-band provides a slightly more strict constraint of the IG dust than V-band, 
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Perlmutter et al. 



we reg ard the upper limit of 0.1 mag as that in the 5-band. In addition 
(1999) report that the difference of the observed reddening between the local and distant 
SNe is (E{B - V)) z=0 . 5 - (E{B - V)) z=om = 0.002 ± 0.03 mag. Although there seems to be 
no systematic difference, we can still consider the absolute value of the colour excess by the 
IG dust less than 0.03 mag. In summary, the IG extinction and reddening from observations 
of distant SNe la at z ~ 0.5 are Af*f =05 0.1 mag and \E(B - V)^ 5 \ < 0.03 mag. 

Suppose an observer who observes a source with a redshift of z at a wavelength of A Q b s 
in his/her rest frame. The amount of the IG extinction is given by 

mag 

where a is the grain radius, Q is the extinction efficiency factor, c is the light speed, and 
n iGU j g jq g ra j n number density in a unit physical volume, which is 

Xfi b p c ,o(l + zfZ{z) 

Ud {Z) ~ W^3 ' (5j 



,086^fg( a ,^)nr(/)( IT ^ M , (4) 



where g(= 2 gem 3 ) is the grain material density. For simplicity, we have assumed that dust 



grains distribute uniformly with corresponding to the redshift. We do not consider any 



structure of the dust distribution. 

In figure 3, we show the amount of the IG extinction divided by \ for a source with 
z = 0.5 as a function of the assumed grain size. While the results for the high SFH is 
depicted in the figure, the extinction amount for the low SFH case is only a factor of about 
1.5 smaller than that of the high SFH case. The extinction amount is independent of the 
grain size as long as 2ira < A. This is caused by (1) the extinction cross section (= Qext^o 2 ) 
is proportional to a 3 because the extinction efficiency factor, Q ex t, is proportional to the 
grain size, a, linearly, and (2) the number density of the grain has a dependence of a -3 for 
a fixed dust mass. On the other hand, Q ext becomes constant (almost 2) when 2ira > A, 
so that the extinction cross section is determined by mainly the geometrical one which is 
proportional to a 2 . Thus, the extinction amount is proportional to a -1 . When 2na ~ A, 
grains interact with photons the most effectively, so that the amount of extinction shows a 
peak in both panels of figure 3. 

We shall remember the observational constraints of A l g^L 05 ^ 0.1 mag and \E(B — 
V r )^2o^ 5 | ^ 0.03 mag. Hence, we can obtain the upper bound of x y i a Eq.(4) or figure 3, 
which is shown in figure 4. The silicate and graphite cases are shown in the panel (a) and 
(b), respectively. The dotted and dashed curves indicate the upper bounds of x based on 
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Figure 3. Intergalactic extinction in various bands against a source with the redshift z = 0.5 as a function of the IG grain 
size. The vertical axis is divided by the parameter x> which is defined by equation (3). The panels (a) and (b) are the cases 
of silicate and graphite, respectively. The star formation history assumed is the lower case shown in figure 1. The solid curves 
indicate the extinction amount in the U, B, V, R, I, J, H, and K bands from top to bottom. 



A l ^^L 05 = 0.1 mag and \E(B — V) l ^^\ = 0.03 mag, respectively. The discontinuity of the 
dashed curve is due to E(B — V) = 0; the left of the discontinuity is positive E(B — V) and 
the right is negative one. The solid curve is the upper bound based on the thermal history 
of the IGM which is obtained in the next section. The rejected area of x is shaded; the 
thin shade means the rejected area based on the SNe la extinction/reddening, and the thick 
shade is the area based on the IGM thermal history (section 4). We show only the high SFH 
case. 

We find that, for both grain types, the upper bound of x from the distant SNe la is 
~ 0.1(a/0.1 jum) for 0.1 /jm < a < 1 /im. While we have no constraint of x f°r the small 
(a < 100 A) silicate grain (panel [a]), the upper bound of x is ~ 0.1 when a < 100 A for the 
graphite case (panel [b]). This difference is caused by the different optical properties between 
graphite and silicate; small silicate is more transparent than graphite in the optical bands. 
For 100 A < a < 0.1 /xm, the upper bound of x shows a local minimum corresponding to the 
peak shown in figure 3. Finally, there is not any constraint of x f° r a ver Y large (a > 1 /xm) 
grain. 
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Figure 4. Rejected area of x> the mass ratio of the IG dust to the total metal as a function of the IG grain size: (a) silicate 
and (b) graphite case. The dotted and dashed curves are the upper bounds from extinction and reddening of distant SNe la. 
The solid curves are the upper bound obtained from the thermal history. Two spectral indices of the background radiation are 
considered; a = 1 (bottom solid) and 2 (top solid). The thick and thin shaded areas are the rejected areas based on the thermal 
history (a = 1 case) and on observations of SNe la, respectively. The high SFH is assumed. 



4 CONSTRAINT FROM IGM THERMAL HISTORY 



As shown in the previous paper ([inoue fc Kamaval 12003), we show that the amount of the 
IG dust is constrained by using the thermal history of the IGM. When a dust grain is 
hit by a photon with an energy larger than a critical value, an electron escapes from the 
grain; photoelectric effect. Such a photoelectron contributes to the gas heating if its energy 
is larger than the mean kinetic energy of gas particles. As shown in Appendix A (see also 



Nath. Sethi, fc Shchekinov 



1999: 



Inoue fc Kamava 



2003), the photoelectric heating by dust 



grains is comparable with, and sometimes dominate, the atomic photoionization heating 
in the IGM. Of course, the efficiency of the dust heating depends on the dust amount. If 
a model of the IGM has too much dust, the theoretical temperature of IGM exceeds the 
observational one owing to the dust photoelectric heating. Therefore, we can put an upper 
bound of the amount of the IG dust so as to keep the consistency between the theoretical 
IGM temperature and the observed one. In the next subsection, we describe how to calculate 
the thermal history of IGM affected by the dust photoelectric heating. An upper bound of 
X, which represents the amount of the IG dust, is estimated in subsection 4.2. 
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4.1 Thermal history of IGM 

In this paper, we consider only a mean t emperature of the IG M, Ttqm, for simplicity. The 



Tjqm time-evolution is described by (e.g. 



Hui fc Gnedin 



19971 ) 



T lGM dX 2(T-A) 

— Z_Hi IG M 77 77" + 



(6) 



dt ' X dt ' ^ksXrib 

where H is the Hubble constant, is the cosmic mean number density of baryon, V and A 
are the total heating and cooling rates per unit volume, respectively, X is the number ratio of 
the total gaseous particles to the baryon particles, i.e., X = J2 rii/n^, where rii is the number 
density of the i-th gaseous species and we consider H I, H II, He I, He II, He ill, and electron. 
We neglect the effect of helium and metal production by stars on the chemical abundance 
for simplicity. Fortunately, time evolution of their abundance is not important. Indeed, the 
metal mass fraction reaches at the most 0.002 (figure 2). The number ratio of helium to 
hydrogen is always about 0.1 after the Big-bang. A constant mass fraction (Y = 0.24) of 
helium is assumed throughout our calculation. 

We solve equation (6) coupled with non-equilibrium rate equations for these gaseous 
species by the fourth-order Runge-Kutta scheme. These rate equations with rate coefficients 
and the heating/cooling rates are summarised in Appendix B. In the calculation, the time- 
step is adjusted to being 1/1000 of —X m / (dXm/dt) for dX m /dt < and Xn U /(dX m /dt) 
for dXfti/dt > (see equation [B6] for dXm/dt). The number density of the IG dust 
at each redshift is determined by equation (1), (2), and (5) depending on \. The grain 
ch arge and heating/cooling rates are determined by a standard manner which is described 



in 



Inoue &: Kamava 



(2003) (see also Appendix A). 



The initial condition is as follows: the s tarting redshift is z 



that the Hell reionization occurred (e.g. 



Theuns et al 



3.4, at which it is considered 



2002ah . The initial temperature is 



set to be 25,000 K, which is the m ean IGM temperature at this redshift suggested by 



the Lyman a forest in QSO spectra (jSchave et al 



2000; 



Theuns et al 



2002H ). We assume 



an ionization equilibrium balanced between the recombination and the photoionization as 
the initial chemical abundance. In each time step, the calculated chemical abundance at 
z < 3.4 reaches almost another ionization equilibrium balanced among the recombination, 
the photoionization, and the collisional ionization. 3 



3 For a technical reason, we did not include the collisional ionization in the calculation of the initial condition; to avoid being 
dXm/dt = for the time-step. Since the collisional ionization plays only a minor role, the calculated chemical abundance is 
different slightly from the abundance in the recombination-photoionization equilibrium. 
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Figure 5. Redshift evolution of ionizing background intensity. Data points with error-bars are taken from lScott et all 120021) . 
The solid line is an analytical fit described in equation (7). 



The background radiation is required to calculate the IGM thermal history. We assume 
a single power-law background radiation; its mean intensity at a frequency of v is J v = 
Jl(z//z/l) _q , where Jl and z/l are the mean intensity and the frequency at the hydrogen 
Lyman limit. We also assume the spectral index a to be constant, but the intensity Jl 
evol ves along the redsh ift. In figure 5, such an evolution is displayed. The data are taken 
from Scott et al. ( 20021 ) who investigate the QSO proximity effect on the number density of 
the Lyman a forest in spectra of QSOs at z = 0-5, and estimate the Lyman limit intensity 
of the background radiation in the redshift range. Here, we use a fitting formula as 

1 \ « 2.5 x 10- 23 (1 + zf* (7) 

erg s 1 cm z sr 1 Hz 1 

for < z < 4, which is shown in figure 5 as the solid line. The background radiation at z < 3 
is likely to be dominated by QSOs. Hence, we consider two cases of a = 1 and 2. Such values 



of a are consistent with the QSO dominated background radiation (jHaardt fc Madau 



Zheng et al 



1996 



1997). In appendix A2, we show that the results with a = 1 single power-law 



background s pectrum should be consi stent quantitatively with those with a more realistic 



spectrum like lHaardt fc Madaul ((1996). 

In figure 6, we show some examples of the IGM thermal history, assuming 0.1 (xm IG 
dust and the spectral index of a = 1. The panel (a) shows the silicate and the high SFH 
case, and the panel (b) shows the graphite and the low SFH case. In each panel, six cases of 
X are depicted as the solid curves. By definition of x ( ec i-[3])> X — means no IG dust and 
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Figure 6. Examples of the IGM thermal history: (a) silicate for high SFH and (b) graphite for low SFH cases. The grain size 
is assumed to be 0.1 um. Six solid curves in each panel are the thermal histories ass umed a value of \, t he mass ratio of the 
IG dust to the total metal, indicated on each curve. The data points are taken from lSchave et all l2000l) . The spectral index 
of the background radiation is assumed to be unity. 



X = 1 mea ns that all of metal is condensed into the IG dust. The observational data are 
taken from Schave et al.l (j200ol ). They observe the Lya clouds with the column density of 



10 13-15 cm -2 (i.e., slightly over density regions), and convert the temperature of the clouds 
into that at the mean density of the IGM by using the equation of state of the IGM. Thus, 
we can compare both thermal histories directly. In the next subsection, such a comparison 
is presented quantitatively. 



4.2 Constraint for \ from thermal history of IGM 



Once theoretical histories of IGM temperature are obtained, the amount of IG dust is con- 
strained from the comparison of the theoretical temperature with observational one. Hence, 
we compare our theoretical thermal histories with 10 observational points at the range of 
1.5 < z < 3.4 of ISchave et al.l ()2000l ). We reject a case of too much \ by the least squares 
method. The rejection criterion is the significance level less than 30% in x 2 -test. The ob- 
tained upper bound of \ in the high SFH case is shown in figure 4 as the solid curves. In 
this figure, the results from distant SNe la are overlaid as the dotted and dashed curves. 
Moreover, the rejection areas from the thermal history and distant SNe la are distinguished 
by the thick and thin shadings, respectively. 

Combining constraints from the thermal history with those from distant SNe la, we 
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Table 1. Upper bounds of x 
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a By thermal history. 
^ By SNe la observations 
c By definition 



obtain the rejected range of x as a function of the IG grain size, which are summarised in 
Table 1. We find a rough upper bound of x as 0.1 with a factor of a few uncertainty, except 
for a very large (~ 1 \xm) case and a medium-size (~ 100 A) silicate of a — 2. 

The upper bound of x from the thermal history has a positive dependence of grain size. 
This corresponds to the fact that the dust heating rate has a negative dependence of grain 
size (see figure A3). Especially, for small silicate grain, we obtain a strict upper bound of 
X from the thermal history, whereas observations of distant SNe la cannot provide any 
constraints. 

We notice here that the upper bound of x f° r small silicate is smaller than that for 
small graphite. This is why a small (a < 0.1|xm) silicate has a larger efficiency factor for 
absorption than that of graphite in the UV-X-ray regime. Hence the grain potential, mean 
photoelectron energy, and heating rate of small silicate are larger than those of graphite 
(figures A1-A3). Moreover, we find a positive dependence of the spectral index a against 
the obtained upper bounds of x, which is due to the negative dependence of a against the 
dust heating rate (figure A3). 
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Figure 7. Maximum (a) IG dust density and (b) dust-to-gas mass ratio in the IGM as a function of rcdshift. In each panel, 
the solid and dashed curves are the cases of high and low star formation histories, respectively. The assumed \, the mass ratio 
of the IG dust to the total metal in the universe, is 0.1. 

5 DISCUSSION 

We have obtained upper bound of x as a function of grain size from observations of distant 
SNe la and comparison of theoretical IGM thermal history with observational one. Here we 
discuss what our results imply. 

5.1 Allowable amount of IG dust 

How much dust can exist in the IGM? Once assuming a value of x, we obtain the IG dust 
density via equation (5). In figure 7, we show the upper bounds of the IG dust mass density 
and V IGM as a function of redshift. The solid and dashed curves are the cases of high and 
low SFHs, respectively. The assumed upper bound of x is 0.1 for both cases of SFHs. The 
uncertainty of this value of x is about a factor of a few as long as the IG grain size is 
smaller than 1 |j.m and the background spectral index a — 1 (see also Table 1). We find 
that the upper bound of the local (z ~ 0) universe is determined well; the local IG dust 
density is less than about 1CT 34 g cm" 3 , or equivalently the dust-to-gas ratio is less than 
about 3 x 10~ 4 which is about 1/20 of the Galactic value. Along the redshift, the allowed 
dust density increases and the dust-to-gas ratio decreases. The increasing/decreasing rates 
change at z ~ 1 at which the (I+2) dependence of the assumed SFH changes (eq.[l]). Taking 
into account a factor of 2 uncertainty of the upper bound of x, we obtain the maximum IG 
dust density as 
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§ cm 3 [ (i±2) (for ^ > 1) 

or the maximum dust-to-gas mass ratio as 



l+z 



-0-7 

(for z < 1) 



« w = (i - 4) x i(r 4 <; ^^^-(1.5-3) ;: . (9) 



2 



(for z > 1] 



hedust-to-gas ratio at z ~ 3 is consistent with the previous result by 
( 20031 ) . 



Inoue fc Kamava 



5.2 IG extinction and reddening 

In figure 8, we show the upper bounds of IG extinction and reddening expected from the 
upper bounds of x m the case of the high SFH and a = 1 summarised in Table 1. Four 
cases of grain size as 10 A, 100 A, 0.1 (xm, and 1 \xm are indicated by dotted, dot-dashed, 
solid, and dashed curves, respectively. The constraints from SNe la observations of Ab ^0.1 
mag and \E(B — V)\ ^ 0.03 mag at z = 0.5 are also shown as the downward arrow in each 
panel. We note that the vertical axes of panels (c) and (d) are the absolute value of colour 
excess. Indeed, colour excesses of 1 (xm for silicate, and of 0.1 \xm and 1 [im for graphite are 
negative. 

We find that the upper bound of the IG extinction is ~ 0.2 mag for a source at z = 1 
from pa nels (a) and (b) of figure 8 . This value agrees well with the result from SDSS quasars 
data bv iMortsell fc Goobarl (J2003). For z > 1 objects, the upper bound of the IG extinction 
becomes ~ 0.5 mag, and as an extreme case, we cannot reject the possibility of 1 mag IG 
extinction for a source at z ~ 3. 

Interestingly, we can investigate the nature of the IG dust by using the IG reddening. 
For z > 1 sources, the expected absolute value of colour excess by the IG grain larger than 
~ 0.1 [im is very small, at the most ~ 0.05 mag, whereas that by a smaller grain can 
reach 0.1 mag or more. Thus, it may be possible to determine a typical size of the IG grain 
from observations of colour excess against a source at z > 1; the detection of > 0.1 mag 
colour excess for such a source proves the existence of small (< 100 A) IG grains. If the IG 
dust is dominated by such small grains, the composition of the IG dust can be found. The 
small graphite grains show a prominent absorption feature at 2175 A. Thus, we expect a 
local minimum of colour excess for a source at z ~ 2.5 as shown in panel (d) of figure 8. 
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Figure 8. Maximum intergalactic extinction and reddening against a source at a redshift in the case of high star formation 
history; observer's S-band extinction of (a) silicate case and (b) graphite case, and colour excesses between B and V-bands 
in observer's rest frame of (c) silicate case and (d) graphite case. Assumed \ values are summarised in Table 1. For dotted, 
dot-dashed, solid, and dashed curves, we assumed the grain size of 10 A, 100 A, 0.1 um, and 1 um, respectively. The spectral 
index of the background radiation is assumed to be unity. The upper bounds from observations of distant SNe la are shown 
as the downward arrow in each panel. The vertical axes of panels (c) and (d) are the absolute value of the colour excesses. 
Actually, colour excesses of 1 um for silicate and of 0.1 and 1 um for graphite are negative. 



Therefore, if we detect such a change of colour excess along the redshift, we can conclude 
that many small graphite grains exist in the IGM. 

Observations to detect the IG extinction and reddening are very challenging but strongly 
encouragin g. High-z gamma-ray b ursts can be good background light sources for such ob- 



servations ([Perna fc Aguirre 



200(1 . 
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5.3 Ejection efficiency of dust from galaxies 

The dust-to-metal mass ratio in the Galaxy, 5mw, is 0.3-0.5 (e.g. 



Spitzer 



1978; 



Whittet 



20031 1. We may consider the dust-to-metal ratio in the IGM, 5 IG m is equal to #mw> ^ ^mw is 



typical for all galaxies, and metal and dust are ejected together from galaxies to the IGM 
keeping <5mw- Is this compatible with the obtained upper bound of x ^ 0.1? 

A fraction of metal produced by stars in galaxies exists out of galaxies. This met al escape 
fracti on is defined as fz, esc - While fz, esc is still uncertain, an estimate of it is 50-75% (jAguirre 



1999L and references therein, see also our figure 2). We shall define another parameter of the 
IGM metallicity; Zjq M . Because of % — T> IGU /Z and /z,esc = Zi GM /Z, where Z is the total 
cosmic metallicity, we find <5igm = ^ igm /^igm = x//z,esc- Unless /z,esc is less than ~ 0.5, 
5\gm is estimated to be smaller than 0.2 if x < 0.1. Therefore, our result of x ^ 0.1 with 
/z,esc ~ 0.5 may indicate that 5 IGM < #mw- 

If that is true, we have to consider some mechanisms to reduce <5igm during the dust 



trans fer. For example, dust destruction during the transfer from galaxies to the IGM ( Ag uirre 



1999), and/or different ejection efficiencies between metal and dust . Time evolu tion of the 



Tuouel fl2 



the dust- 



dust-to-metal ratio in galaxies may be also important. As shown by 
to- metal ratio in younger galaxies (i.e., higher- z galaxies) may be much smaller (~ 70% off) 
than the present value of the Galaxy. In the case, a time-averaged <5mw can become smaller 
than the current 5mw adopted above, so that our constraint of x ma J t» e cleared. In any 
case, we cannot obtain a rigid quantitative conclusion at the moment, because uncertainties 
are still large. Further studies of this issue are very interesting. 



5.4 IGM temperature at low redshift 

As shown in Nath et al. (1999) and Appendix A, the dust photoelectric heating becomes 
more efficient for a lower gas density. Although the background intensity decreases along 
the redshift (figure 5), the decrement of gas density is more efficient than the decrement 
of the background intensity, so that the importance of the dust heating increases for a 
lower redshift. While we obtained constraints of the amount of the IG dust from the IGM 
temperature at z ~ 2-3 in section 4, the IGM temperature at a lower redshift of z < 
1 provides us with a further constraint of the IG dust. Therefore, to measure the IGM 
temperature at z < 1 is very interesting. 

Here, we demonstrate how much temperature is allowed by our upper bounds of x- Figure 
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Figure 9. Maximum IGM temperatures corresponding to the upper bounds of \ in table 1 for the case of silicate, background 
spectral index of unity, and high SFH. The d otted, dot-dashed, a nd solid curves are the case of no IG dust, 100 A, and 0.1 um, 
respectively. The data points are taken from lSchave et all J2000l) . 

9 shows the IGM thermal histories assumed the upper bounds of \ m Table 1 for the case of 
background spectral index a = 1 and high SFH. The dotted, dot-dashed, and solid curves 
are the case of no IG dust, 100 A silicate, and 0.1 (xm silicate, respectively. We can make a 
very similar figure for graphite case. The temperatures shown in the figure are upper bounds, 
which is denoted as Ttgm,u P - 

After checking all cases listed in table 1, we find that for a smaller (< 100 A) grain case, 
except for graphite of a = 2, Ttgm,u P at z < 1 is still much higher than 10,000 K. On the 
other hand, for a larger (> 0.1 yim) case, except for 1 yim graphite with high SFH and a = 1, 
^igm,up at z < 1 becomes lower than 10,000 K as well as no IG dust case. Therefore, we may 
conclude that IG grains are small if temperature higher than 10,000 K is observed at z < 1. 
Conversely, a lower IGM temperature at z < 1 provides us with a very strict constraint 
against small IG dust. 

6 CONCLUSION 

We investigate the amount of the IG dust allowed by current observations of distant SNe 
la and temperature of the IGM. The allowed amount of the IG dust is described as the 
upper bound of \i the mass ratio of the IG dust to the total metal mass in the Universe. 
To specify Xi two models of cosmic history of metal production rate are assumed. That is, 
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we have assumed two cosmic star formation histories expected from the recent observation 
of the high redshift objects. Our conclusions £1X6 clS follows: 

(1) Combining constraints from the IGM thermal history with those from distant SNe 
la observations, we obtain the upper bounds of x as a function of grain size in the IGM; 
roughly x < 0.1 for 10 A < a < 0.1 fj.m, and x < 0.1(a/0.1 \im) for 0.1 \im < a < 1 p.m. 

(2) The upper bound of x ~ 0.1 corresponds to the upper bound of the IG dust density; 
the density increases from ~ 10~ 34 g cm~ 3 at z = to ~ 10~ 33 g cm~ 3 at z ~ 1, and keeps 
a constant value or slowly increases toward higher redshift. 

(3) The expected IG extinction against a source at z ~ 1 is less than ~ 0.2 mag at the 
observer's _B-band. For higher redshift sources, we cannot reject the possibility of 1 mag 
extinction by the IG dust at the observer's 5-band. 

(4) Observations of colour excess against a source at z > 1 provides us with information 
useful to constrain the nature of the IG dust. If we detect ~ 0.1 mag colour excess between 
the observer's B and V-bands, a typical size of the IG dust is < 100 A. Moreover, if there 
are many graphite grains of a < 100 A in the IGM, we find a local minimum of the colour 
excess of a source at z ~ 2.5 corresponding to 2175 A absorption feature. 

(5) If half of metal produced in galaxies exists in the IGM, the obtained upper bound of 
X ~ 0.1 means that the dust-to-metal ratio in the IGM is smaller than the current Galactic 
value. It suggests that some mechanisms to reduce the dust-to-metal ratio in the IGM 
are required. For example, dust destruction in transfer from galaxies to the IGM, selective 
transport of metal from galaxies, and time evolution of the dust-to-metal ratio in galaxies 
(i.e., a smaller value for younger galaxies). 

(6) Although we obtain constrains of the IG dust from the IGM temperature at z ~ 2-3, 
the temperature at z < 1 provides us with a more strict constraint of the IG dust. For 
example, the detection of temperature higher than 10,000 K at z < 1 suggests that the IG 
dust is dominated by small (< 100 A) grains. 
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APPENDIX A: DUST PHOTOELECTRIC HEATING IN IGM 

The dust photoelectric effect in the IGM is summarised. The basic equ ations of the dust 
photoelectric effect are presented in various places, for example, section 2 in llnoue fc Kamava 
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(2003). More detailed information on this effect can be found in IWeingartner fc Draine 



mm - 



We consider spherical silicate and graphite grains. Under a condition suitable for the IGM 
(low density and intense UV radiation), these grains have a positive electric charge, which 
is determined mainly by the competition between the collisional electron capture and the 
photoelectric ionization. The ion collision is negligible, while the proton collision is included 
in the calculation in section 4 and for making the following figures. The charge on grains is 
in an equilibrium state, which is achieved quickly (~ 10-100 yr). 

The input parameters to obtain the equilibrium charge are grain type, grain size, gas den- 
sity, gas temperature, radiation intensity, and radiation spectrum. We assume the incident 
radiation spectrum to be a power-law. In figure Al, we show these dependences of the equilib- 
rium grain potential energy normalized by the gas kinetic energy, i.e., eU jk^T = Z^e 2 fak^T, 
where U is the grain potential, is the grain charge, a is the grain size, T is the gas tem- 
perature, e is the electron charge, and &b is the Boltzmann constant. We examine three 
cases of the power-law index of the incident radiation, a = 1, 2, and 5, where the mean 
intensity is J v and proportional to v~ a . While the spectrum of the incident radiation is 



likelv to be a power-law with index of 1-' 


Haardt & Madau 


1996 




Zheng et al. 


1997) 



mated by QSOs (e.g., 

reference of a very soft background radiation. The parameter set assumed in the calculation 
are noted in each panel of figure Al; a t _5 = a/0.1(xm, n t _5 = n/10~ 5 cm -3 , T 4 = T/10 4 K, 
and J -2i = J UL /10~ 21 ergs -1 cm -2 sr -1 Hz -1 , where z/l is the Lyman limit frequency. These 
values may be suitable for the IGM at z ~ 3. 

The dotted lines in figure Al (a) and (c) indicate an upper limit of the grain potentia l 
based on an estimate of the critical potential for the ion field emission (jDraine fc Haol l2002): 
eU jk-oF < 3500(T/10 4 K) _1 (a/0.1 (j.m). If the grain potential exceeds this upper limit, singly 
charged ions may escape one by one from the grain surface, so that the grain is destroyed 
gradually. For panels (b) and (d), this upper limit is out of the panels. We can conclude that 
this process is not so important for our interest. 

In figure A2, we show a mean photoelectron energy from dust grains normalized by the 
gas kinetic energy, (E ve ) /k^T, as a function of (a) grain size, (b) gas density, (c) gas temper- 
ature, and (d) radiation intensity. Moreover, the ratio of the dust photoelectric heating rate 
to the hydrogen photoionization heating rate (r po /TH P i) is depicted as a function of these 
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quantities in figure A3. For the dust heating rate, we assume the dust-to-gas mass ratio to 
be 10 -4 as a nominal value, which is about 1/100 of the Galacti c dust-to-gas mass ratio. 



Nath. Sethi, Shchekinov 



The obtained results are roughly consistent with those of 
()1999l ) . While some quantitative differences are seen between our results and theirs, they 



may be caused by differences of the adopted photoelectric yield, absorption efficiency factors, 
and radiation spectrum. 

In summary, we find that the dust photoelectric heating exceeds the photoionization 
heating for a case of a smaller grain size, a lower gas density, a higher temperature, a more 
intense radiation, or a harder radiation spectrum. Furthermore, we find that silicate and 
graphite grains show similar results, and all of eUjk^T, (E pe )/kBT, and r pe /TH P i show a 
power-law like dependence of grain size, gas density, gas temperature, and radiation inten- 
sity. These behaviors will be derived analytically with some approximations in the next 
subsection. 



Al Analytic investigation 

To understand the behavior of grain charge and other quantities found in figures A1-A3 ana- 
lytically, we will adopt s ome further approximatio ns in this subsection whic h do not appear in 



the method described in 



Inoue fc Kamava 



Inoue fc Kamava 



( 2003T ). However, the method by 
(2003) is used in the calculation for making figures A1-A4 and in section 4. We would like 
to ask the reader to be careful in this point. 

Let us express eU jh^T as x. If we neglect the charging rate by grain-ion collision, the 



charge equilibrium is expressed as 

ax 



s e n e (v e )(l + x) 



Q V Y V — — dv 
hv 



(Al) 



where s e is the sticking coefficient for an electron collision, n e is the electron number density, 
(v e ) is the mean kinetic velocity of electron, Q u is the absorption efficiency factor of grains, 
Y v is the photoelectric yield, and J v is the mean intensity of the incident radiation. The 
left-hand side of the above equation is the electron capture rate per unit area, and the right- 
hand side is the photoelectric ionization rate per unit area. Although the integral should 
be summed up to infinity in general, we set the upper limit of the incident photon energy, 
lii/ mH = 1.24 keV. We also adopt s e = 1. 

In order to inves t igate analytically we adopt a simple functional form for the photoelec- 
tric yield bv iDraind (|1978r ): 
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Figure Al. Normalized equilibrium potential of grains as a function of (a) grain size, (b) gas number density, (c) gas tempera- 
ture, and (d) mean intensity at the Lyman limit. The solid and dashed curves indicate silicate and graphite cases, respectively. 
For the incident radiation, we consider three cases of the power-law index, J„ oc v~ a \ a = 1, 2, and 5 from top to bottom 
curves in each panel. Assumed values for parameters are indicated in each panel; a,_s = a/0.1nm, 71,-5 = n/W~ 5 cm -3 , 
T, 4 = T/10 4 K, and J_ 2 i = J„ L /1CT 21 erg s 1 cm 2 sr 1 Hz 1 , where v\, is the Ly man limit frequency. The dotted lines in 
panels (a) and (c) indicate an upper limit of the grain potential iDraine fc H ao 2002). 



Y = Y 



hp 



(A2) 



for hv > I p and Y v — for otherwise, where I p is the ionization potential; I p ~ W + xk-gT 
with W being the work function . We no te that a m ore realistic functi o n of t he photoelectric 



yield by 



Weingartner fc Draind (j200 1b) is used in 



Inoue fc Kamava 



( 2003T 1 and in section 



4. Moreover, we adopt an approximation form of the absorption efficiency factor as Q u ~ 
Qhi^/^h)' 13 , where z/l is the Lyman limit frequency of hydrogen. For a grain larger than 
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Figure A2. Mean energy of dust photoelectron normalized by gas kinetic energy as a function of (a) grain size, (b) gas number 
density, (c) gas temperature, and (d) mean intensity at the Lyman limit. Notations are the same as figure Al. 



~ 0.1 pm, (3 — 0, and for a smaller grain, (3 ~ 1-2 against ultraviolet photons. A power-law 
spectrum for the radiation, J u = J\Xy / Vl) _c * is also assumed. 

If we define a function as f(hu) = QYJjhv, df/dhu = only when hv = hv* = 
I p (a + (3 + 2) /(a + (3 + 1). Because the function / has the single peak at hv*, the integral 
in equation (Al) is an order of (4ir/h)Q(hv*)Y(hi / *)J(hi'*). If x ^> 1, then, equation (Al) 
is reduced to 

l+a+ p 47rQ L y oo J L f a + (3+l \ a+P fhv h \ Q+/3 /^rm^X V* 

~ s c n c h(a + (3 + 2) \a + (3 + 2) \k B T ) \8k B T J ' 1 ' 
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Figure A3. Ratio of dust photoelectric heating rate to hydrogen photoionization heating rate as a function of (a) grain size, 
(b) gas number density, (c) gas temperature, and (d) mean intensity at the Lyman limit. Notations are the same as figure Al. 
Additional parameter, T> is the dust-to-gas mass ratio; £>,-4 = X>/10~ 4 . 



where we approximated J p ~ xk-gT because Wjk^T ~ 1 for T ~ 10 4 K, and substituted 
(v e ) = {%kv,T /ixm c ) l l 2 . If n c ~ n with n being the gas number density, we find 
TT / T n ( \ \ 1 /( 1 + a +P) 

eU_ ^ / JiUl[0) \ r _ ( i /2+Q+/ 3)/(l + a+/3) _ ^ 

k-gT y n J 

Indeed, such dependences are found in figure Al. Moreover, Q L ~ (a/O.lpm) for a < 0.1 
\xm and Ql ~ 1 otherwise. Thus, eU jk^T shows nearly no dependence of grain size for a 
large size. 

The mean photoelectron energy is defined as 
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_ / E pc Q u Y u AttJ u / hvdv 
K pc) fQuYA-n Jul 'hvdv ' 1 } 

where E pe is the energy of the photoelectron. We can express E pe « i](hu — 7 P ), where rj is a 

numerical factor less than unity beca use a part of energy of the inc ident photon is converted 



into the phonon energy of the grain ((Weingartner &: Draine 



2001b). 



f we adopt a parabolic 



functio n for the energy distribution function of the photoelectron as IWeingartner fc Draine 
(200rbJ), the numerical factor rj = 1/3-1/2 depending on the energy of the incident photon. 
We adopt rj — 1/2 below. 

If we define a function as g{hv) = E pc QYJ/hv and assume the functional forms adopted 
above for Q, Y, and J, dg/dhv is zero only when hv = hv** = I p (a+/3+2)/(a+P). Thus, the 
integral in the numerator of equation (A5) is an order of E pe (hv**)Q(hv**)Y(hv**) J(hv**). 
If x ^> 1, we obtain 

I • (A6) 



(E pe ) 2 / a + P 



xk B T a + (3 \a + (3 + 1 J 
Therefore, (E pe )/k-QT has a similar parameter dependence to ell/k^T, which is observed in 
figure A2. 

The saturation of (E pe ) /k^T is seen in the case of a = 1 for a lower density in figure A2 
(b). This may be due to the effect of the maximum photon energy assumed in the calculation. 
We do not consider the incident photon energy higher than 1.2 keV. For the case, I p reaches 
several hundreds eV, so that the peak energy hv** becomes nearly the maximum energy. 
Thus, the above estimate may be an overestimate for such a case. 

The total photoelectric heating by dust is given by r pc = n<j7, where is the grain 
number density and 7 is the heating rate per a grain. While rid oc na~ 3 for a certain dust-to- 
gas ratio, 7 is proportional to (-E pe )a 2 nT 1//2 a; because the photoelectric ionization rate per a 
grain balances with the electron capture rate per a grain, where a 2 dependence comes from 
the geometrical cross section of grains. We have assumed x ^> 1 again. On the other hand, 
the photoionization heating rate, TH P i, is proportional to n 2 T~ ' 7 in the photoionization 
equilibrium, where the temperature dependence comes from the recombination coefficient. 
Therefore, we find 

^oca-VT 2 - 2 , (A7) 

iHpi 

where we have used the relation (E pc ) oc xT (equation [A6]). Remembering J, n, and T 
dependences in x described in equation (A4), we can understand J, n, and T dependences 
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shown in figures A3 (b), (c), and (d). Because Ql oc a for a < 0.1 |am and Ql ~ 1 for 
otherwise, we see a double power-law dependence of a in panel (a) of figure A3. 



A2 Effect of the spectral break 

As shown by lHaardt &: Madaul (J1996), the real spectrum may show a break at the He II 
Lyman limit (54.4 eV), whereas we have assumed a spectrum without break. Here this point 
is discussed. Assuming the power-law spectral index (a) is fixed all over the spectral range 
for simplicity, we multiply the intensity of the background radiation above the He II Lyman 
limit by a factor of /ken, which is called the spectral break factor in this appendix. As 
extreme cases, fueii = 1 means that there is no break, and /neii = means that there is no 
photon above the He II Lyman limit. In all calculations, except for figure A4, /heii = 1 has 
been assumed. 

Figure A4 shows the effect of fueii on the normalized mean photoelectron energy from 
dust grains {(E pc ) / 'k^T) , which indicates the heating efficiency per a grain. Only the silicate 
case is shown, but the graphite case is very similar. For a = 1 case (dotted curves), we 
observe the photoelectron energy (E ve ) /k^T decreases with decreasing /neii- This is because 
the number of high energ y photons decrease s if the spectral break is larger (i.e. smaller fneii)- 



According to figure 5 in 



Haardt fc Madad (|l996h . the break is significant like fueii ~ 0.1. 



One might think that our assumption of fueii = 1 with a = 1 results in an overestimation 
of the dust heating, so t hat a larger amount of d ust may be allowed in the IGM. 

However, figure 5 of Haardt fc Madau ( 19961 ) also shows that the spectral index is 0.5 
rather than unity for < 1 keV photons which we are interested in. In the a = 0.5 case 
(solid curves of figure A4), we find a good quantitative agreement with the case of fneii = 1 
and a = 1 (top dotted curve) if fueii = 0.1-0.3. Therefore, our results obtained from the 
background spectrum of /neii = 1 and a = 1 should be quantitatively very consistent with 
those from a more realistic spectrum with the He II Lyman limit break. 



APPENDIX B: CHEMICAL RATE EQUATIONS AND COEFFICIENTS 

In this paper, we consider HI, HII, Hel, Hell, Helll, and electron as gaseous species. That is, 
we neglect the effect of the metal production by stars. The primordial helium mass fraction 
Y = 0.24 is always adopted throughout our calculation. 

Let us define a non-dimensional number abundance of each gaseous species as 
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Figure A4. Effect of the spectral break at the He II Lyman limit (54.4 eV) for silicate grains. Solid and dotted curves are the 
spectral index a = 0.5 and 1 cases, respectively. For each type of curves, we assume the spectral break factor, /h ii = 0.1, 0.3, 
and 1 (no break) from bottom to top, respectively. The gas density n__5 = n/10 -5 cm -3 , the gas temperature T4 = T/10 4 K, 
and the intensity of the incident radiation J -21 = J^/IO -21 ergs -1 cm -2 sr _1 Hz -1 are assumed. 



V - Ui 
A-i = — , 

n b 



(Bl) 



where n« is the number density of i-th species, and n D is the baryon number density, which is 
given by rib(z) = f2bPc,o(l + z) 3 / /x with p Cj0 being the local critical density, and \i being the 
mean mass of baryon particles. Baryon is assumed to be only hydrogen and helium. Using 
the helium mass fraction Y, hence, we obtain 
4(1 -Y) 



Y 



4-3F ' 



X m + X mi - 4 _ 3y , 

-^Hel + -^Hell + -^Helll = 

and 

^ = v4^w) mp ' 

where m p is the proton mass. For electron abundance, we have 



X e — -Xffll + -^Hell + 2X He III 



Chemical rate equations for the gaseous species are 
dX m 



dt 
rfX Hc i 
dt 

and 



= X H iiX e n b aHii - X H i(X c n b /9Hi + 7hi) , 

= -^Hell^e^b^Hell — -^Hel (^e^bAlel + 7Hel) , 



(B2) 
(B3) 

(B4) 

(B5) 

(B6) 
(B7) 



32 A. K. Inoue & H. Kamaya 



Table Bl. Recomb ination and collisional ionization coefficients in cm 3 s 1 taken from ICenl Il992f) and 
iTheuns et all Il998i) . 



Recombination 
ami = 6.30 x 10- 11 T- 1 /2 r -0- 2 /(i + t 6 7 ) 

QHdi = 1-50 x 10- 10 T-°- 6353 +ag eII 
ttHein = 3.36 x 10- 10 T- 1 / 2 T 3 - V(1 + T 6 ' 7 ) 
Dielectric recombination 
"Soil = L9 x 10- 3 T- 1 - 5 e- 4 ' 7 x 1 ° 5 /T(l + o.3e- 9 ' 4 x lo4 /T) 
Collisional ionization 
Pm = 5.85 X 10- 11 T 1 /2(l + T l/2)-l e -157809.1/T 
/3 H eI = 2.38 X 10~ 11 T 1 /2(1 + T l/2)-l e -285335.4/T 

feell = 5.68 x 10- 12 T 1 /2(1 + r 5 1/2 )- 1 e-63isi5/T 



Table B2. Parameters for ionization cross sections from lOsterbrockl Jl98 



species 


ff< [10~ 18 cm 2 ] 


i>i [10 15 Hz] 


h 


Si 


HI 


6.30 


3.29 


1.34 


2.99 


Hel 


7.83 


5.94 


1.66 


2.05 


Hell 


1.58 


13.2 


1.34 


2.99 



dXn in 

j£ = — ^Helll^e^baHelll + -^Hell (-^e^bAlell + 7HeIl) , (B8) 

where a*, and 7« are recombination coefficients, collisional ionization coefficients, and 
photoionization rates for i-th species, respectively. The adopted functions of a« and /3j are 
tabulated in Table Bl. The photoionization rate is given by 

H = / a i)U ——diy « -7-<TiJ L — — — — , B9 

where J u = JiX v l v \)j a with z/l is the hydrogen Lyman limit frequency, and a^ v is the 
ionization cross section: 



for v ^ Vi and otherwise er^ = (lOsterbrock 



(BIO) 



1989). The parameters for ionization cross 



sections are summarised in Table B2. The last term in equation (B9) is valid when z/ max ^> j/j. 
In this paper, we set hv mBX = 1.24 keV. 

By solving equations (B6)-(B8), we obtain Xm, Xn e j, and Xhgiii- Once these frac- 
tional abundances are obtained, other abundances are found from equations (B2), (B3), 
and (B5). In addition, the numbers of hydrogen and helium nuclei are constant since the 
helium mass fraction Y is now constant. Thus, dXmi/dt = —dXm/dt and dXn e n/dt = 
—dXn e i/dt — dXu e iii/dt. Therefore, the term of dX/dt in equation (6) of section 4.1 is re- 
duced to — dXm/dt — dXn c i/dt + dXn c m/dt because of dX/dt = J2idXi/dt = dX e /dt and 
equation (B5). 

We now consider the atomic photoionization heating and the photoelectric heating by 
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Table B3. Coo l ing c oefficients in erg cm 3 s 1 taken from ICenl 119921) . iTheuns et alj il998T) . and 
llnoue fc Kamaval <2003tl . 



Recombination cooling 
HII 
Hell 
HcIII 
dust a 

Dielectric recombination cooling 
Hell 

CoUisional ionization cooling 
HI 
Hel 
Hell 

CoUisional excitation cooling 
HI 
Hell 

Bremsstrahlung cooling 
ion'' 

Inverse compton cooling 
CMB photon c 



8.70 x 10- 27 T 1 / 2 T 3 "°- 2 /(l + T°- 7 )X c X ml 

1.55 x 10- 26 T°- 3647 Js: o XhoII 
3.48 x 10- 26 T 1 / 2 T 3 - ' 2 /(l + T 6 7 )X c X Ho iii 
(3/2)k B TR c n d /nl 

1.24 x 10 -13 T -i.S e -4.7xlo B /T( 1 + o.3e- 9 - 4xl ° 4 /T) XcXHc 

1.27 x 10- 21 T 1 /2(i + T 5 1/2 )- 1 e- 157809 - 1 /T XeXm 

9.38 X 10- 22 TV2(1 +T l/2 ) -l e -285335.4/T XcXHci 

4.95 x 10- 22 TV2(1 +T 5 1/2 )- 1 e- 631515 /Tx e X H oii 

7.5 x 10- 19 (1 + T 5 1/2 )- 1 e- 118348 / T X c X H i 
5.54 x 10- 17 T-°- 397 (l + T 5 1/2 )- 1 e- 47638 / T X c X Hc ii 

1.42 x 10- 27 9f T 1 /2x o (X H „ + XhoII + 4X HcIII ) 

5.406 x 10" 36 [T - T CMB (1 + «)](! + z) 4 X a /n h 



a Re is the electron capture rate per a grain and is the grain number density. These quantities are calcu- 
lated by the method described in llnoue fc Kamaval J2003ft . Although the factor (3/ 2) is an approximation , 



this electron capture cooling is not important. A more realistic case is investigated in lDraine fc SutTnl Il98' 
" gi is the Gaunt factor and assumed to be 1.5. 

c Tcmb is the current temperature of the cosmic microwave background and set to be 2.7 K. 

dust as the heating mechanism of gas in equation (6). The photoionization heating rate per 
a z-th species atom/ion is given by 



[hv-hvijOi^— — dv « 47rcTi^J L — 
hv V vt 



2h - 1 



+ 



h - 1 



,(bh; 



ya + Sj — 1 a + Si a + Si + 1^ 
where parameters, <7j, Ui, bi, and Sj are summarised in Table B2, and a is the power-law 
spectral index of the incident radiation. Again, the last term of equation (Bll) is valid when 
^max v i- The dust photoelectric heating is given by calculating t he equilibrium cha r ge an d 



the ejection rate of the photoelectron by the manner described in llnoue fc Kamava (2003). 
Finally, the adopted cooling rates are summarised in Table B3. The metallic line cooling is 
not so important for our problem because the tempe rature interested is less tha n 25,000 K 



and the metallicity is 1/10-1/100 of the Solar value ((Sutherland fc Dopital ll993). 



